Phenomenology of Jet Quenching in Heavy Ion Collisions 
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We derive an analytical expression for the quenching factor in the strong quenching limit where 
the pt spectrum of hard partons is dominated by surface emission. We explore the phenomenological 
consequences of different scaling laws for the energy loss and calculate the additional suppression of 
the away- side jet. 



m 
o 
o 



m 

> 
oo 
m 
o 

00 
O 
<N 
O 



o 

=5 



X 



It is commonly believed that the yield for "hard" ob- 
servables in high-energy nuclear reactions scales as the 
number of binary nucleon-nuclcon (NN) collisions occur- 
ring during the encounter of the two nuclei. This expec- 
tation applies to processes that are characterized by a 
high virtuality q 2 , for which final state interactions are 
negligible, such as lepton pair production or the total 
yield of heavy flavor quarks. There is no a priory reason 
to expect this rule to hold for high-g 2 processes, in which 
the final state can be strongly modified by interactions 
with comovers. A well documented example is the pro- 
duction of heavy vector mesons, such as the J r /\E f , which 
is found to be "anomalously" suppressed in collisions of 
heavy nuclei at the CERN-SPS Q. 

The yield of hadrons produced with high transverse 
momentum px in Au+Au collisions at the Relativistic 
Heavy Ion Collider (RHIC) has recently been shown to 
be significantly suppressed in comparison with the cumu- 
lative yield of NN collisions |2,|3|- This effect, called "jet 
quenching" , was predicted to occur as a result of energy 
loss by the hard scattered partons due to interactions 
with the surrounding dense medium 0,0,0. The theory 
of this energy loss has been a topic of intense research 
over the past few years 0, B IS lii! 111! ■ The present con- 
sensus is that the dominant mechanism for the energy 
loss in QCD is collisionally induced radiation of gluons 
by the fast parton. 

It is difficult to measure the energy loss of a scattered 
parton directly in heavy ion reactions, because the large 
multiplicity of emitted hadrons makes it almost impossi- 
ble to isolate the resulting jet by kinematic cuts. How- 
ever, the energy loss of the parton is imprinted as an 
equivalent loss of en ergy of the leading hadron produced 
in its fragmentation 12]. This is what has been observed 
in the RHIC experiments. Generally, it is assumed that 
the fragmentation occurs after the parton has left the 
comoving medium, and thus is described by the mea- 
sured vacuum fragmentation functions. We will, there- 
fore, not be concerned with the conversion from partons 
to hadrons, but focus directly on the p^-spectrum of scat- 
tered partons. 

Preliminary data from Run 2 at RHIC (Au+Au at a 
center-of-mass energy of 200 GeV per nucleon pair) con- 
firm the effect observed in Run 1 and its interpretation 
as jet quenching 0. For pions with pt ~ 5 GeV/c the 
measured suppression factor Q{pt) is about 1/5. There 



is not yet complete agreement between different exper- 
iments about the scaling of the hadron yields with the 
number of nucleons participating in the reaction, iV part , 
or with the binary collision number N co \y. The PHOBOS 
collaboration has presented evidence that the yield of 
charged hadrons withp^ ~ 4 GeV/c scales like N valt 
This finding is contradicted by data from the PHENIX 
collaboration [J^. It is not clear whether the discrep- 
ancy is due to different normalization methods, different 
experimental acceptance, or other reasons. 

In addition, the experiments have found that the su- 
pression factors for mesons and baryons are quite differ- 
ent up to transverse momenta of 5 Gev/c. It was re- 
cently proposed that this phenomenon can be attributed 
to a competition between different hadron formation pro- 
cesses |l6l Il7| , with parton recombination dominating at 
lower pt and fragmentation at higher pr- The partici- 
pant or binary collision scaling may be complicated in the 
transition region. The situation is predicted to simplify 
for pt > 6 GeV/c, where hadrons are overwhelmingly 
produced by fragmentation of an energetic parton. 

We will first show analytically that the spectrum of 
high-pT partons, in the limit of large energy loss, is dom- 
inated by partons emitted from the surface of the colli- 
sion zone and thus scales like the surface rather than the 
volume of the interaction region (lRq . We will then ex- 
plore the scaling of the quenching factor Q(pt) with par- 
ticipant number and Pt- We finally calculate the addi- 
tional suppression of the away-side jet and the azimuthal 
anisotropy of the parton yield in noncentral collisions. 

We begin by considering the loss of energy by an ener- 
getic parton traversing a homogeneous, static medium of 
thickness L. We assume that the geometry is given by a 
cylinder with radius R, as in the boost-invariant Bjorken 
model figj j for a nuclear collision with impact parameter 
6 = 0, and that the parton moves in the transverse plane 
in the local rest frame of the medium. We further assume 
that the effective energy loss, defined as the shift of the 
momentum spectrum of fast partons, depends on pt and 
L in the following general way: 



Ap T = r\p T L . 



(1) 



where /i is a scaling exponent. The linear dependence on 
L holds when the loss occurs in subsequent, independent 
interactions with the medium. It has also been shown 
to be valid, when multiple interactions are suppressed by 
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the LPM effect, due to the steep fall-off of the parton 
spectrum with px fiof . The traversed pathlength inside 
the medium for a parton created at transverse position r 
with an angle <p relative to the radial direction is 



L(cf>) = (R 2 



r 2 sin 2 i 



a/2 



2i?C0S< 



(2) 



where r = Irl, z 2 = R 2 



and the approximation 



is valid near the surface (z <C R). Perturbative QCD 
predicts that the parton s pec trum at moderatedly large 
values of px has the form 
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with a power v w 8 and po 
spectrum is given by 



1.75 GeV/c. The quenched 
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(4) 



Replacing the integration over r with an integration over 
z, using the approximation (J2J, and formally extending 
the range of the integration to infinity, one finds: 



Q(pt) 



2(po + Pt) 



ixRr]{y 



(5) 



Two things are remarkable about this result. First, the 
factor R in the denominator reduces the scaling of the 
parton yield with the size of the reaction zone by one 
power of i?, from a volume to a surface dependence. Sec- 
ond, the dependence of Q{pt) on px is determined by 
the power fj, governing the pr-dependence of the energy 
loss. For /j, = 1 the quenching factor Q falls slowly with 
increasing px] for \l = 1/2 it grows with px, implying 
less quenching at higher px ■ 

We have confirmed the range of validity of the approx- 
imate analytical expression JSJ by comparing it with the 
exact integral Q. As seen in Fig. the analytical ap- 
proximation Qanai deviates from the exact result by less 
than 5% when the quenching factor Q < 0.2. An even 
better agreement is found when the analytical result is 
divided by the correction factor (1+Q 2 nal ). The excellent 
agreement suggests that we may extend the calculation to 
noncentral collisions. Generalizing the surface-to- volume 
ratio of the cylindrical geometry (S/V = 2/R) to the ge- 
ometry of a collision with impact parameter b, one finds 



Q(PT,b) = Q(p T ,0)a b /(a b - sina b ) 



(6) 



with a b = 2 arccos(fe/2_R). This approximation is valid 
with about the same accuracy as (O, except for very 
peripheral collisions. 

In order to be able to address the experimental data, 
we need to relax some of the geometrical oversimplifica- 
tions used above. 
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FIG. 1: Comparison between the analytical approximation 
© and the exact result @ for the quenching factor Q. The 
curves are for Au+Au at b = and pr = 6 GeV/c, with 
a — 1/2. The calculation assumes a homogeneous transverse 
profile for jet production and quenching medium. Also shown 
is the expression Q a nai/(1 + 1 



?anal ) 



1. The transverse profile of the primary jet yield is 
not homogeneous, but proportional to the binary 
NN collision profile 



T(r,b) = Pl (r)p 2 (r-b) 



(7) 



where pi (r) is the longitudinally integrated density 
of nucleus i. 

2. The density of the comoving medium is also not 
homogeneous in the transverse plane. We assume 
that it is proportional to the local density of par- 
ticipant nucleons which, in the Glauber model, is 
given by 

Ppart(r) = Pl (r)(l- e -^ r - b )) 

+p 2 (r-b)(l- e - CT " 1 W) , (8) 
where a denotes the inelastic NN cross section. 



3. The comover medium expands and its density de- 
creases with time. Here we assume 



p(r,r) = Cp part (r)/(r + r ) 



(9) 



which is modeled to represent a longitudinal, boost 
invariant expansion. As an approximation to the 
phenomenology at full RHIC energy we use the val- 
ues C « 3 and tq — 1 fm/c. 

Finally, we need to return to equation for the en- 
ergy loss. Perturbative QCD predicts that the radiative 
energy loss depends quadratically on the medium thick- 
ness L As pointed out by Baier et al. 01 > this holds 
for the average energy loss AE of a given parton, but 
the average energy loss of observed partons with fixed 
transverse momentum px has a different scaling. This is 
so, because the energy loss distribution D(e) is strongly 
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skewed toward small values of e by the steeply falling px- 
spectrum of fast partons. In fact, the average shift of the 
spectrum due to the en ergy loss, here called the effective 
energy loss, is given by jlfj| 

Ap T « ot s L\J irqpx / v , (10) 

where q encodes the "scattering power" of the medium, 
which is proportional to the density. For an expanding 
medium, the expression qL 2 must be replaced with 

q L eS = 2q J o rdr—^, (11) 

where r(r) = r + vr denotes the position of the fast 
parton in the medium at time t, and go is a function of 
the transverse position r at which the jet is produced. 
We thus can make contact with by writing 

Ap T = r]'L e &*Jp(r,0)p T /v (12) 

with the constant rf — a s ^7rq / p(0), which does not 
depend on r. 

We will denote the scaling law 112JI for the energy loss 
as BDMS. In our following numerical study we have ex- 
plored two other scaling laws. The first one is the Bcthc- 
Heitler (BH) scaling law [UH^ 

Ap T = i]p T dr p(v(t),t) ee r]p T (Lp) e g (13) 
Jo 

corresponding to )i = 1. The second scaling law is 

Ap T = r\p T \/ {Lp) cS , (14) 

which we will call the (RW) scaling law. It could be 
interpreted as describing a random walk in px as the 
fast parton traverses the medium, with some interactions 
resulting in an energy gain and others in a loss of energy. 

We begin the discussion of our numerical results for 
the quenching factor Q with its dependence on the trans- 
verse momentum of the fast parton, shown in Fig. EI The 
QCD-motivatcd BDMS law (solid line) and the other two 
scaling laws exhibit clearly different behaviors. This re- 
flects the different pt scaling of the energy loss in these 
models (linear for BH and RW; square-root for BDMS). 
The data from the RHIC experiments 0, 0, E| suggest 
that the quenching first becomes stronger with increas- 
ing momentum, reaches a minimum, and finally begins 
to diminish. This would indicate that the BDMS law 
only applies at high pt, and that other laws govern the 
energy loss at lower px [HI or hadron production is not 
dominated by parton fragmentation in this kinematic re- 
gion. We note that the dependence of Q on px is well 
described by the analytical formula JHJ). 

The impact parameter dependence of the quenching 
factor is shown in Fig. [31 plotted as the yield per half 
the number of particpant nucleons against the particpant 
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FIG. 2: Dependence of the quenching factor Q on px for 
central collisions. The parameter r\ is chosen such that 
Q(Pt) ~ 0.2 for pr = 10 GeV/c in each case. The scaling 
laws (BH, RW) exhibit stronger quenching with increasing 
Pt, in agreement with preliminary RHIC data, in contrast to 
the BDMS law. Equation J^J provides a good description of 
the dependence on px seen here. 
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FIG. 3: Quenched hard parton yield divided by half the num- 
ber of participant nucleons as a function of iVpart for Pt = 10 
GeV/c. The values of the stopping power strength parameters 
are r? = 0.06 (RW), 77 = 0.017 (BH), and 77' = 1.1 (BDMS). 



number N part . The unquenched jet yield, scaling with the 
number of binary NN collisions, would increase relative 
to -ZVpart. As the figure shows, the quenching counteracts 
this increase, and the yield per participant actually falls 
for the BDMS and the BH laws as the collision centrality 
increases. An approximately flat behavior, as observed 
in the PHOBOS experiment ^ij, is only found for the 
RW scaling law. 

For noncentral collisions, the quenching factor Q is a 
function of the azimuthal emission angle, because the ge- 
ometry is not axially symmetric. This is known to lead 
to an angular asymmetry of a quadrupole shape in the 
spectra of high-p-r particles [23|,|24(. The elliptic flow pa- 
rameter i>2 is defined as the Fourier component propor- 
tional to cos(20) of the angular distribution of particles 
with respect to the scattering plane [i^. We find (see 
Fig. 0} that the values of v-i < 0.1 obtained for all three 
scaling laws are significantly smaller than the measured 
values (i>2 ~ 0.2) for semicentral and peripheral colli- 
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FIG. 4: Elliptic flow parameter V2 as a function of the number 
of participants for the three energy loss models. The values 
of the parameters as the same as in Fig. [3] 
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FIG. 5: Incremental away-side jet suppression factor Q aa j as 
function of the "same-side" jet suppression factor Q for the 
BDMS energy loss law. 

sions [2(| . However, the calculated V2 for partons would 
be large enough to explain the measured elliptic flow of 
hadrons, if the hadrons were produced by recombination 
in the p^-range of the RHIC data • 

As observed hadrons from hard partons preferentially 
originate from the surface region facing the detector, the 
parton emitted in the opposite direction has to traverse 
more material and thus endures an even larger energy 
loss. This leads to an additional suppression of the away- 
side jet and its leading hadron spectrum. The depen- 
dence of the incremental away-side hadron suppression 
factor <2 as j on the primary suppression factor Q is shown 
in Fig. [3] for the BDMS scaling law for two impact pa- 
rameters (6 = and b = 8 fm). There appears to be a 
universal relationship, which is linear for Q > 0.2. 

Several conclusions can be drawn from our results. 
First, the momentum dependence of the BDMS energy 
loss formula does not seem to agree with some of the 
RHIC data. A linear dependence of the average energy 
loss on pt is in better agreement with the data in the pr 
region explored so far [21| . Also, the dependence of the 
effective momentum loss Apx on the medium thickness 
L e g predicted by BDMS does not yield a scaling of the 
charged hadron yield with participant number as seen 



in the PHOBOS data up to p T = 4.25 GeV/c; only the 
RW scaling law yields such a dependence. This may in- 
dicate that energy gain and loss mechanisms are compet- 
ing in this momentum range |2sj . Another explanation 
may be that hadrons at intermediate momenta are not 
produced by fragmentation of fast partons, but by other 
processes, such as parton recombination. The observed 
magnitude of the elliptic flow lends support to this in- 
terpretation. Finally, we have found a universal relation 
between the same-side and away-side suppression factors 
for the BDMS law, which can be tested experimentally. 
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